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Abstract 

It is shown that highly localized solitons can be created in lower dimen- 
sional Bose-Einstein condensates (BEC), trapped in a regular harmonic 
trap, by temporally varying the trap frequency. A BEC confined in such 
a trap can be effectively used to construct a pulsed atomic laser emitting 
coherent atomic wave packets. It is also shown that one has complete 
control over the spatio-temporal dynamics of the solitons. The dynamics 
of these solitons are compared with those constructed in a BEC where the 
trap frequency is constant. 

Pacs(03.75.-b, 05.45.Yv, 03.75.Pp) 

1 Introduction 

The creation of the Bose-Einstein condensate (BEC) [1] has opened up an area 
with a potential for huge technological application. On one side, active research, 
both theoretical and experimental, is being done to understand the properties 
of this new phase of matter and on the other side practical applications of the 
condensate are being developed. The fact that a BEC constitutes coherent 
matter waves analogous to the coherent light waves, makes it the most suitable 
candidate for the construction of an atom laser. The first atomic laser was con- 
structed in 1997 in MIT [2] and other experimental results were reported in [3J, 
[4 . An atom laser can be used to develop accurate atomic clocks [5] , gyroscopes 
for ultra precision navigation [5] and will also help in the development of co- 
herent atom optics [7], [5J, atom interferometry [S], [TU] and atomic transport 
[11] . The realization of bright solitons and soliton trains [12] -[18] in a BEC has 
given further impetus to this line of research. 

Solitons are ubiquitous in nature and have a unique property that they do 
not disperse even after traveling long distances. Optical solitons have been 
successfully used for the purpose of telecommunications [TU]. A bright soliton 
is a bunch of atoms moving together and its density is more compared to the 
background. It has been shown in [20 that solitonic pulses can be generated in a 
BEC in an elongated harmonic trap, by simultaneously changing the scattering 
length and by making the trap expulsive. Generation of soliton pulses for a 
pulsed atom laser, using the quantum state transfer property was discussed by 
Liu et. ai, [21] . The dynamics of an atom laser were studied in [22] - [24] 
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It is a well known fact that the Gross-Petaivskii (GP) equation describes the 
dynamics of a BEC at the mean field level. An analytical method to solve the 
GP equation for a BEC trapped in an elongated harmonic trap was presented in 
|25j . Here, both the scattering length and the trap frequency were time depen- 
dent and the phase of the solitonic solution got related to the trap frequency 
in a non-trivial way. This relation was exploited to obtain complete control 
over the spatio-temporal dynamics of the condensate, in addition to obtaining 
exact expressions for the solutions. This became possible because this relation 
provided one with the freedom to vary the trap frequency and also to tune the 
atomic s-wave scattering length. This in turn allowed one to study the effects of 
variations in trap frequency on the dynamics of the BEC and also to see how the 
nonlinearity changes with time. In [26], the effect of sudden temporal changes 
in the trap frequency of an expulsive harmonic trap have been investigated and 
highly localized solitons were observed in certain cases. 

In the present paper, the emphasis is on the variations of the trap frequency 
of a regular harmonic trap and their effect on the soliton dynamics. It is shown 
that in a BEC trapped in a regular harmonic trap, suitable temporal variations 
of the trap frequency lead to the creation of highly compressed solitons. The 
trap frequency is varied using an oscillatory potential which asymptotically goes 
to zero. When compared to the solitons obtained with constant trap frequency, 
these solitons are highly localized. From the context of a pulsed atomic laser, 
highly localized solitons are desirable. Hence, the concentration will be on the 
generation of localized bright solitons and it will be shown that these solitons 
are stable for certain ranges of t. Therefore, one studies the evolution of the 
bright solitons in these specific regions only. 



2 Solutions of the NLSE 

The dynamics of a dilute Bose gas trapped in a cylindrical harmonic trap with 
a time dependent confinement in the z direction given by V = m'to^x 2 + 
y 2 )/2 + m'w 2 (i)z 2 /2, is governed by the GP equation. This equation is reduced 
to the quasi-one dimensional nonlinear Schrodinger equation (NLSE), with a 
tight confinement in the z direction (lu± » uj z ), using a Guassian trial wave 
function [55], [27], [25], 

% d t i> = ~d„i/> + 7(*)M 2 V> + ^M(t)*V, (1) 

where j(t) = 2a s (t)/aB, M(t) = uj^(t) / u>\ and a± — h/m'uj±. Here, 7(t) is 
the nonlinear coupling which is attractive if j(t) < and repulsive if 7(i) > 0. 
a s is the scattering length and as is the Bohr's radius. Note that M (t) > 
corresponds to the regular oscillator and M(t) < corresponds to the expulsive 
oscillator. The exact solutions of the NLSE are obtained by using the ansatz 
solution 

iP(z,t) = y/AtyF{A(t)[z-l{t)]}exp[t*(z,t)], (2) 
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where 

<$>(z,t)=a(t)-\c{t)z 2 (3) 

with a(t) = ciq + L A 2 (t')dt', clq being a constant. Here, A(t) gives the 
amplitude and l(t) gives the location of the center of mass of the soliton respec- 
tively. Substitution of ip(z, t) in the NLSE gives the elliptic differential equation 
for F 

F"(T)-\F{T) + 2kF 3 (T) = (4) 

with k = T = A(t)(z — l{t)) and the differentiation is with respect to 

T. Solutions of this equation are the Jacobi elliptic functions [35], namely 
sn (T, m), cn (T, to) etc with elliptic modulus m. Importantly, solitonic solutions 
arise in the limit m — ► 1. Specifically, in the attractive nonlinear regime (70 < 
0), F(T) = cn(T/ro,TO) in ((2|) and it describes a bright soliton train, with 
Tq = — Aom/^fo, A = (1 — 2to)/Vq. In the limit to — > 1, cn (T, to) — > sech(T) 
and the same equation describes a bright soliton 



i>(z,t) = v/I(Osech(T/To)expHt) - -c(t)z\ (5) 

with Tq = — Aq/jq. Similarly, the dark soliton trains are obtained with F(T) = 
sn(T/ro,m) in ([2]) in the repulsive nonlinear regime (70 > 0) with Tq = Aq/jq, 
A = — (1+to)/tq . The dark soliton is obtained in the limit to — ► 1 as sn (T, to) — ► 
tanh(T) and is given by 



ip(z, t) = y/A(i) tanh(T/r ) exp[ia(t) - -c(t)z 2 }. (6) 

Substitution of © in the NLSE flTJ| gives us the following consistency conditions, 
with A and 70 being constants, 

A(t) = A a exp (J* c(t')dA , 7(0 = 7o^, (7) 

and 

^-+c(t)l(t)=Q. (8) 

which in turn give us expressions for the control parameters. It is also found 
that c(t) satisfies the Riccati equation 



dc(t) , 



dt 



c z (t)=M(t), (9) 



which can be mapped to a second order differential equation using the change 
of variable 

c(t)=-|h#(i)]. (10) 

The new equation, 

- 0"(t) - M(t)ct>(t) = 0, (11) 
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is in the form of the Schrodinger equation, if one writes the ratio of the trap 
frequencies as M(t) = Mq + V(t). Here, V(t) is like a time dependent potential 
with eigenvalue Mo- This allows one to modify the temporal profile of the trap 
frequency using potential models for which the Schrodinger equation is exactly 
solvable. Thus, if one knows the solutions <p(t) of (fTTj) , one can obtain c(t). 
This when substituted in (O and fl5|), gives the consistency conditions in terms 
of 4>{t) as follows 

A(t) ^„||, 7(t) . 70 |), i(t) . io ||, (12) 

where Iq is a constant. Thus one obtains exact expressions for the control pa- 
rameters. These allow one to write the complete solution of the NLSE using 
which one can study and control the spatio-temporal dynamics of the solitons. 
Using (|9|) and (|10p one can vary the trap frequency and obtain corresponding 
soliton solutions and study their spatio-temporal dynamics. 



3 Soliton dynamics 



As mentioned in the introduction, the present paper looks at the soliton dy- 
namics when the trap frequency of the regular harmonic trap is modelled by 
an oscillatory function which vanishes as t — > 0. Such a variation in the trap 
frequency is obtained by taking the deformed free particle potentials [30] as 
V(t) in (fTTj) . The eigenfunctions of these potentials are known explicitly and 
are used to obtain the solitonic solutions. 
Case I The first potential to be used is V(t) = V(t, A), t > 0, 



V(t,X) 



32M cos 4 (VMot) 8M sin(2^/M 't) 



n 2 



Do 



(13) 



where D = 2^/Mo~t + sm (2v^Mo i) +4\/Mo A, VAio > and A > 0. The square 
integrable solution for pip with the above V(t) is 



& (*, A) 



4v r Mocos(v r Moi) 
Do 



with energy eigenvalue Mq. Using the above solution in (p~2|) . one obtains 

D 



and 



A(t) = A 
lit) = 7o 



4AVMo 

4A % /A/o 
4A/Mo" 



D 



sec(VM i), 
sec(v/Moi) 

COS (\J Mq t). 



(14) 

(15) 
(16) 

(17) 
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Substitution of the above equations in ([5]) gives the bright soliton solution. 

The variation of the nonlinearity with time and the dynamics of the solution 
as it evolves in space and time are shown in the fig. 1. It can be seen from 
fig. 1(a) that the nonlinearity changes sign periodically and the change is dis- 
continuous. Since the emphasis is on the study of bright solitons, one chooses a 
range where the nonlinearity is negative and continuous as shown in fig. 1(b). 
In this range of t, the matter wave density \ip(z 7 t)\ 2 , is evolved for varying z 
and t as shown in fig. 1(c). One can see that the bright soliton is localized and 
with time gets compressed. The amplitude of the soliton can be controlled by 
varying Aq in (|15[) and the location of the soliton can be varied by changing Iq 
in (TTT)) . Thus, one gets to control the dynamics of the soliton completely. For 
comparison we plot the dynamics of the bright soliton when the trap frequency 
is constant i.e. V(t) — 0. Thus (fTTj) is nothing but the Schrodinger equation for 
a free particle. We use the solution <fii(t) = cos (/Mo t) and obtain the bright 
soliton solution and plot its evolution which is shown in fig. 1(d). 

When one compares the dynamics of these two solitons taking the same 
parameter values, one observes that in the earlier case the soliton is highly 
localized and compressed whereas in the later case, the soliton moves within the 
specific t range with more or less the same amplitude, undergoing little or no 
change. Also one observes huge amplification of the soliton in the case where 
the trap frequency is time dependent. It has been see that for both the cases, 
smaller the value of y/Mo, bigger will be the range of t in which the nonlinearity 
is continuous, being either attractive or repulsive and hence longer the life time 
of the soliton. From the above example it is clear that by varying the trap 
frequency temporally, one can generate solitons which are highly localized and 
amplified. Hence, this kind of traps can be used in the construction of pulsed 
atom lasers. 

In the above case, the solitons have been evolved for t > 0. One can do 
a similar analysis for t < with A < and it turns out there is no change 
in the dynamics of the soliton. The only difference lies in the behavior of the 
nonlinearity plotted in fig. 2(a). One can see that within a certain range of time 
and with the same parameter values, the nonlinearity changes sign smoothly 
as shown in fig. 2(b). This implies that within this range the nonlinearity 
changes its character from being attractive to being repulsive in a smooth way 
without one resorting to the use of the Feshbach resonance [31] - [34] . This 
when compared to the case where the trap frequency is constant, one observes 
that the regions of attractive and repulsive nonlinearity get interchanged in the 
two cases and j(t) changes its nature periodically and discontinuously along the 
negative t axis for the later case. 

Case II Another example for which the soliton dynamics is studied is when the 
trap frequency is modulated using 

.... 32Af sin 4 (VAioi) 8M sin 2 (2 VMq t) 

V 2 (t,X) = - x (18) 

for t > 0. Here, D x = - sin (2 -/Mo t) + 4-/MoA, /Mo > and A > 0. 

This potential is also an oscillating potential which vanishes at infinity. The 
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square integrable solution of (fTTj) with the above V(t) is 



hit, A) = 



4Vlio sin [y/Mot) 



(19) 



with eigenvalue Mq > 0. Thus using the above solution in (|12p . one obtains 



A(t) = A Q 



D\ sin (y 



D[ sin (a 



and 



7(<) = 



L>i sin (-/Mo*) 
3 Z^sin (VMoM' 



(20) 
(21) 

(22) 



where D[ = 2/Mo *o — sin (2/Mo *o) + 4/MqA. As in the previous case, the 
nonlinearity changes sign periodically and discontinuously. 

Hence, for a certain range of t where j(t) is continuous and negative, the 
soliton is evolved and the dynamics are shown in fig. 3(a). For comparison the 
dynamics of a soliton, for the same parameter values, in a BEC located inside a 
trap with constant trap frequency i.e., (V(t) = 0) are given in figure 3(b). Here 
4>2{t) = sin (y/Mo t) was used to obtain the analytical expression for the bright 
soliton. In contrast to case I, both the solitons evolve in the same way and have 
almost the same amplitude except that in the case where the frequency is time 
dependent, the soliton is more localized. 



4 Conclusions 

In the present paper, it has been shown that in a BEC confined is a regular har- 
monic trap, a suitable temporal variation of the trap frequency, allows one to 
construct highly localized and compressed solitons in the scenario where there is 
no loss or gain of atoms in the condensate. The difference in the soliton dynam- 
ics when the trap frequency is varied and when it is constant can be seen very 
clearly in case I. Thus, by using suitable oscillatory, asymptotically vanishing 
potentials to vary the trap frequency, one can construct coherent atomic wave 
packets in a BEC. These can be used effectively in an atomic laser. Moreover 
this type of variation in the trap frequency can be implemented with ease in an 
experimental set up. The fact that the solutions of a linear Schrodinger equation 
can be combined with the solutions of the NLSE effectively to obtain analytical 
solutions of the GP equation allowed one to change the trap frequency. This 
method also gives one a complete control over the soliton dynamics. Finally one 
would like to point out here that one can also construct and study the evolution 
of bright and dark solitonic trains and the dark solitons using the equations 1(5]) 
and ([6]) in appropriate limits. 
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(a) The nonlinearity 7(4) (b) 7(4) in a specific range of t: 5.85 < 

t < 9.36 




(c) Bright soliton in the trap with vary- (d) Bright soliton in the trap with con- 

ing trap frequency stant trap frequency 



Figure 1: The plots 1(a) and 1(b) describe the behaviour of nonlinearity when 
the trap frquency is time dependent (j(t) (y-axis), t (x-axis)). The plots 1(c) 
and 1(d) show the evolution of the bright solitons, in the approximate range 
t — 5.85 to t = 9.36, described by the equation ([5]) (|i/> 2 (z, i)| (z-axis), t (x-axis) 
and z (y-axis)). The parameter values are 70 = —0.5, M Q — 0.83666, A = 1, 
lo = 5, To = 1 and A = 1. 




(a) The nonlinearity 7(4) (b) The nonlinearity 7(t) in a specific 

range —5.2 < t < —1.5 



Figure 2: The plots of the nonlinearity j(t) in the negative t domain {^{i) 
(y-axis), t (x-axis)). The parameter values are 70 = —0.5, Mq = 0.83666 and 
A= 1. 
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(a) Bright soliton in the trap with vary- (b) Bright soliton in the trap with con- 

ing trap frequency stant trap frequency 

Figure 3: The evolution of the bright solitons discussed in case II ( | -0 2 (^;, i) | (z- 
axis), t (x-axis) and z (y-axis)). The range of Ms < t < 3.745. The parameter 
values are 70 = —0.5, M — 0.83666, Aq = 1, Iq = 5, tq = 1 and A = 1. 



10 



